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Quasiperiodic long range order is intermediate between spatial periodicity and disorder, and the 
excitations in ID quasiperiodic systems are believed to be transitional between extended and local- 
ized. These ideas are tested with a numerical analysis of two incommensurate ID elastic chains: 
Frenkel-Kontorova (FK) and Lennard- Jones (LJ). The ground state configurations and the eigenfre- 
quencies and eigenfunctions for harmonic excitations are determined. Aubry's transition by breaking 
the analyticity is observed in the ground state of each model, but the behavior of the excitations 
is qualitatively different. Phonon localization is observed for some modes in the LJ chain on both 
sides of the transition. The localization phenomenon apparently is decoupled from the distribution 
of eigenfrequencies since the spectrum changes from continuous to Cantor-set-like when the inter- 
action parameters are varied to cross the analyticity-breaking transition. The eigenfunctions of the 
FK chain satisfy the "quasi-Bloch" theorem below the transition, but not above it, while only a 
subset of the eigenfunctions of the LJ chain satisfy the theorem. 



PACS numbers: 63.20.Pw; 6f .44.Fw; 6f .44.Br 

There are surprisingly many examples of quasi-one- 
dimensional (ID) quasiperiodic incommensurate sys- 
tems, including quasicrystals, charge density waves, or- 
ganic conductors, and various atomic monolayers ad- 
sorbed on crystalline substrates. Even high-Tc supercon- 
ductors show incommensurate ID-modulation. Although 
embedding these systems in 3D-space may change the 
outcome, understanding the conditions for localization 
of excitations of the underlying model ID-systems is an 
essential beginning. 

Both electron and phonon localization in ID incom- 
mensurate systems have been studied rather extensively 
in recent years. Despite considerable efforts and sub- 
stantial success with several specific models, the prob- 
lem of formulating necessary and sufficient conditions for 
localization is still far from being solved. The underly- 
ing difficulty seems to be related to the specific position 
quasiperiodic order takes among various spatial order- 
ings: it is intermediate between periodicity and disorder 
0. On the one hand, since quasiperiodic systems lack 
periodicity, conventional Bloch theory does not apply. 
On the other hand, the degree of disorder in quasiperi- 
odic systems is not sufficient for Anderson's arguments 
about localization in ID disordered systems ||] to be ap- 
plied to the case. In fact Levitov showed in a rather 
general way that the eigenfunctions of the Schrodinger 
equation for a chain of atoms change from extended to 
localized when the "complexity" of the atomic arrange- 
ments in the chain evolves from periodicity to full-fledged 
disorder. Quasiperiodicity is the boundary between the 
two regimes, and the eigenfunctions in a ID quasiperiodic 
system may be localized, extended, or neither of the two. 
In the latter case the eigenfunctions are called critical. 
As described in the following paragraphs, examples of 
all three behaviors are found in various ID quasiperiodic 
systems. That is, the outcome is model- and system- 
dependent; even minor modifications to model Hamil- 



tonians or variations in parameters may cause drastic 
changes in the eigenfunction behavior. Thus, character- 
izing the solutions of a ID Schrodinger equation with a 
quasiperiodic potential is a difficult problem to solve an- 
alytically. In this paper we take a numerical approach 
to demonstrating the dependence of localization on the 
model and on the parameters. Since there is little differ- 
ence in the formal descriptions of electron and phonon 
localization, we study the latter, because it is easier to 
adapt it to computer calculations without losing essential 
mathematical details. 

An example of extended wave functions in a ID qua- 
siperiodic system was first given in the Dinaburg-Sinai 
theorem [Q for the differential Schrodinger equation 

-^ + Vix)^P = E^P, (1) 

where V{x) is a weak analytic quasiperiodic potential. 
The theorem states that (i) the spectrum is continuous 
and (ii) the wave functions are extended: 

M^) = Mxy, (2) 

where the Wannier functions Uk{x) are quasiperiodic. 
Unfortunately, mathematical constraints, mainly in the 
required analyticity of V{x), do not allow applying the 
theorem to many systems of practical interest. Another 
exactly solvable example is the finite difference (discrete) 
Mathieu equation: 

tjjn+l - 2^n + i^n-l +V{n)lpn ^ Elpn, (3) 

V{n) = V cos(wn -I- a), (4) 

with irrational lo and arbitrary phase a. Aubry and An- 
dre ||] proved that for F < 2 almost all the wave func- 
tions are extended and of the quasi-Bloch form Eq.(||). 
These two examples indicate that sometimes the Bloch 
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theorem may be applicable and the states be extended 
even when the Hamiltonian lacks periodicity. 

The discrete Mathieu equation, Eqs.(||) and (||), also 
gives an example of localization in quasiperiodic systems: 
for y > 2 all wave functions are localized and the 
spectrum is discrete: 



(5) 



Critical behavior, i.e., neither localized nor extended, 
can be found in ID quasicrystals. They are usually mod- 
eled by the discrete Schrodinger equation Eq.(3) with a 
quasiperiodic crystal potential that takes only two values 
representative of two atomic species, A and B, constitut- 
ing the quasicrystal: 

V{n) = V{ujn + a), 
V{x + l) = V{x), 
V[x)=Va, ifO<a:<w 

Vb, ifcj<a;<l. (6) 

This equation was extensively studied by many authors 
1^-^ who gave rather convincing arguments of critical 
behavior; the mathematical proof was obtained by Bel- 
lisard et al. The spectrum is singular continuous, a 
zero measure Cantor set of non-zero fractal dimension. 
The wave functions are neither localized nor extended; 
moreover, they are even divergent: 



mx)VY dx 



\i'{x)\'^ dx 



cx L^, 



(7) 



where 7 is non-integer and eigenenergy-dependent and L 
is the length of the periodically repeated chain cells. 

We believe that solving Eqs.(||) and (^) does not con- 
stitute a general solution for ID quasicrystals. Eq.(||) 
with two undeformable atomic species, A and B, is too 
crude a model even for a ID quasicrystal chain. Since the 
hard sphere approximation leads to the potential V{x) in 
Eq.(^ taking only two values, one for each species, this 
in itself determines the absence of localization for the 
model. It is a theorem that if a quasiperiodic potential 
V{x) takes only finite number of values then localization 
is impossible . To include the possibility of localiza- 
tion, an elastic atomic chain should be examined. 

A simple model embodying many of the features of the 
above systems is the ID Frenkel-Kontorova (FK) chain: 



i(x„+i - Xn- af 



Fcos(27rx„) 



(8) 



Aubry and LeDaeron showed that the minimal en- 
ergy configurations are periodic when the interatomic dis- 
tance a is commensurate with the substrate period and 
quasiperiodic when a is irrational. However, there are 
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FIG. 1. Phase space trajectories for the FK and LJ chains. 
The connected KAM curves (small spots) and discontinuous 
Cantorus curves (large spots) are the trajectories below and 
above the transition by breaking of analyticity, respectively. 
The Cantorus results are shown for varying chain sizes and are 
offset from the KAM curve for clarity. The only discernible 
difference, as a function of the cell length L, is the addition 
of spots at x„ mod 1 = 0.67 and 1.0 from L = 89 to 377 in 
both models, (a) FK chain with V = 0.0207 and V = 0.0271 
for the small and large spots, respectively; (b) LJ chain with 
Vg — —80 and Vg — —110, respectively, and with positions 
normalized to Rg = 3.25 A. 



qualitatively different configurations for the two types 
of incommensurate states separated by the transition by 
breaking of analyticity predicted by Aubry. It is cus- 
tomary to display this transition by plotting trajecto- 
ries in the phase space (a;„_|_i — x„; x„ mod 1) as in 
Fig. 0. For each irrational a there exists a critical value 
Vc of the substrate potential amplitude. U V < Vc, the 
trajectory derived from the atomic configuration {xn} 
is a smooth analytic Kolmogorov-Arnold-Moser (KAM) 
curve; if > the atomic configuration {xn} is dis- 
continuous, and the trajectory is represented by a so- 
called Cantor set torus or Cantorus. The critical value 
Vc depends on a in a rather complicated manner |p^ , |l3| . 
However, it is known that the largest Vc — 0.02461... is 
achieved when a is equal to the Golden Mean |l4| : 



(9) 



We have restricted ourselves to this particular value of a 
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in the numerical solutions. 

The transition, despite its pure mathematical appear- 
ance, has profound physical implications. The smooth 
KAM configurations keep a vestige of translational sym- 
metry: the chain can slide without any change in energy. 
The discontinuous configurations are, on the other hand, 
pinned by the substrate: there is a non-zero Peierls- 
Nabarro barrier |l5|. Thus, the transition by breaking 
of analyticity coincides with an intrinsic pinning transi- 
tion. The spectrum of harmonic excitations also changes 
there: in the KAM regime the lowest eigenfrequency is 
zero (there is a translational Goldstone mode), whereas 
in the pinned Cantorus regime it is non-zero |15| . In view 
of these facts, one is tempted to speculate that some nor- 
mal modes themselves evolve from extended to localized 
when the system is driven over the transition by breaking 
of analyticity. Our numerical results show only a limited 
correlation between the threshold for localization of the 
excitations and the analyticity-breaking transition: the 
correlation does not hold for a chain with anharmonic 
couplings and only a weak localization arises for the har- 
monic FK chain, Eq.(||). 

We checked two chains for phonon localization: 
the Frenkel-Kontorova (FK) and the nearest-neighbor 
Lennard- Jones (LJ). The former is described by Eq.(||), 
the latter by a similar Hamiltonian with only a modifi- 
cation of the interatomic forces ]l^ : 

H = ^[J{xn+i - Xn) + 2K,cos(27rx„/i?s)], 

n 

J{x)^Ae[{alxr-{a/xn (10) 

The results presented here are limited to the LJ chain 
with nearest-neighbor interactions, but tests with the in- 
teraction extended to second and third neighbors showed 
no significant differences. Critical values Vc for the LJ 
chain depend on a as well as on the mean misfit. 

The first step of the solution was to determine min- 
imum energy configurations for Fibonacci approxima- 
tions to the golden-mean misfit of both chains for 
various choices of parameters. For the FK chain, a 
force-relaxation method was used, following guidance by 
Aubry [ p^|Jl5[| on how to avoid getting trapped in lo- 
cal minima: for an initial configuration with equidistant 
atoms, with the first atom in the "proper" place, the sys- 
tem relaxes to the absolute minimum. For the LJ chain, 
both the force relaxation and a gradient search method 
were used to locate the minimum energy configuration 
as there was an additional complication of possible frac- 
turing of the chain 19 1. The transition by breaking 
of analyticity is clearly observable for both the FK and 
LJ chains, as shown in Fig. |l|. We show the smooth 
KAM curve for the largest unit cells, L = 987 and 377 
for the FK and LJ chains, respectively. The clumping 
of the Cantorus is shown there for several values of L 
to demonstrate that its character is well-established at 
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FIG. 2. Integrated density of states for a chain size 
L = 377 for (a) the FK chain and (b) the LJ chain. In (a), 
the integers at the band edges are the count of the eigenen- 
ergies plotted from left to right; the right-most numbers of 
each band are Fibonacci numbers. This coincidence with Fi- 
bonacci numbers also applies to (b). Note the narrow width of 
the bands in (b) for the LJ model as compared with the bands 
in (a) the FK model. Potential parameters and labelings are 
as in Fig. hi. 



modest values of L. 

In the second stage, tridiagonal dynamical matrices 
were constructed for small-amplitude vibrations about 
the ground state configurations {Xn} obtained in the first 
stage. This leads to the following eigenvalue problem for 
the FK chain for unit cell modes of zero wavevector and 
atomic mass m 



V{n) = (27r)Vcos(27ra;„), E = muj^, 



(11) 



and to a similar result for the LJ chain, where the 
eigenenergies E are related to the frequency hy E = 
muPa'^/A^e . Eqs.(pT[) are an example of a general 
quasiperiodic Schrodinger equation. However, unlike the 
case of Eq.(^), the crystal potential V{n) may take an 
infinite number of values. Thus, possible localization is 
not immediately evident pO| ] . The eigenvalues and eigen- 
vectors for the dynamical matrices were obtained using 
standard (EISPAC) bisection routines. 
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FIG. 3. Participation ratio as a function of eigenenergy for the FK chain, (a) full spectral range for L = 377; (b)-(d) over 
a selected range of eigenenergies that corresponds to the central band of (a) for L — 89, 377, and 987. Note the increase in the 
length of the "tails" of the disconnected Cantorus curves as the length of the chain cell increases. Potential parameters and 
labelings are as in Fig. |l|(a). 



The minimal eigenvalue for both the LJ and FK chains 
changes from zero to non-zero when the system is driven 
over the transition by breaking of analyticity in accor- 
dance with earlier observations by Aubry . For both 
models, the eigenvalue distribution also changes over 
the transition: at y < Vc, the eigenvalues are evenly 
spaced, densely filling the bands as in the periodic case; 
at V > Vc, however, the distribution is more fragmented. 
The resulting Integrated Density of States (IDOS) vs. 
eigenenergy E (i.e., squared frequency tu^) and the eigen- 
value distribution over the energy scale are qualitatively 
similar for the FK and LJ chains, as shown in Fig. ||. 
However, the apparent band widths for the LJ case are 
much narrower, and surprising differences arise upon ex- 
amining the spatial character of the eigenfunctions. 

Before presenting the results for localization, we re- 
view some of the limitations of a numerical study of the 
problem. Localization cannot be decided by examining 
the eigenvectors of the dynamical matrix of one rational 
approximation to the Golden Mean misfit. We followed 
the evolution of results for a series of rational approxi- 
mants converging to the Golden Mean. Some properties 
of the eigenfunctions showed strongly convergent behav- 
ior. There were also values of parameters and parts of the 
eigenvalue spectrum where limitations of accuracy and 
slow convergence combined to leave the situation unde- 
cided. This is the case of our rather weak evidence for 
critical behavior in localization. 

Generally, it is difficult to draw quantitative conclu- 
sions about the localization of eigenfunctions by visual 
inspection because of the volume of information involved 
and because oscillatory behaviour is sometimes mixed 
with the decay of a weakly localized state. The data were 



first analysed for evidence of localization by computing 
Participation Ratios (PR)- 



PR=- 



L 



(12) 



If the eigenfunction is extended, PR tends to a finite 
limit as L — !■ cx). If an eigenfunction is truly localized, 
as in Eq.(|^), L x PR tends to a finite limit as the length 
L of the unit cell tends to infinity. Critical states might 
have more bizarre scalings, as in Eq.(|^) po[| . 

We present the results of our localization studies on 
the Frenkcl-Kontorova chain first; participation ratios are 
shown in Fig. ||. In the KAM regime V < Vc, i.e., when 
the ground state is smooth and unpinned, the PR test 
shows that all eigenfunctions are extended. This is cor- 
roborated by demonstrating that the eigenfunctions can 
be transformed to the quasi-Bloch form Eq.(||). More 
precisely, any quasiperiodic function (e.g., Ufc(n)) may 
be represented as: 



Uk{n) = Uk{Lun + a), 
Uk{x) =Uk{x + l). 



(13) 



The periodic generating function Uk {x) in Eq. ( |13[ ) is pre- 
sented in Fig. ^ for three quasimomenta that bracket 
the band shown in Fig ||(c). AX V < Vc, the eigenvalues 
for L ^ oo become doubly degenerate, as in the usual 
application of the Bloch theorem: E{k) — E{—k). We 
determined the value of the quasimomentum k in Eq. (j^) 
by matching to the eigenfunctions of the nearly degen- 
erate eigenvalue pairs. We were able to do this both 
below and above the transition. Below, the function 
is smooth. Above, the eigenvalue pair typically differs in 
the fourth decimal place, but k-values could still be iden- 
tified easily. The function Uk above the transition has a 
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FIG. 4. Quasi-Bloch functions Uk (x) for the FK chain for 
L = 377 for three quasimomenta k and configurations below 
and above the transition, V = 0.0207 and 0.0271, respectively. 
Parameters and labelings are as in Fig. |l|(a). The values of 
k 1^ and the corresponding eigenenergies are (a) k — 1.95, 
E = 1.72 and 1.87; (b) 2.57, 2.13 and 2.22; and (c) 3.11, 
2.51 and 2.55. The quasi-Bloch functions that are shown cor- 
respond to states from (a) the left tail (b) the middle peak 
and (c) the right tail of Fig. ^c). Note the functions above 
the transition maintain the rough shape of the corresponding 
functions below the transition but have a smearing that is 
correlated with the trend of the participation ratios. 
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FIG. 5. Participation ratio for the LJ chain over eigenen- 
ergies representative of the (a) "extended" and (b) "local- 
ized" regimes. Shown are results for L = 55 (diamonds), 
89 (squares), 233 (triangles) and 377 (inverted triangles). 



The open and filled symbols correspond to Vg 



-80 and 



Vg = —110, respectively, for the various chain lengths (these 
are the same potentials as in Fig. 0(b)). The participation 
ratio is shown in (a) whereas the participation ratio times 
the chain length L is shown in (b) to emphasize the local- 
ization over this range of eigenenergies. Note the splitting of 
the curve in (b) as formed by the renormalized participation 
ratios for the various chain lengths. 



clumping along the abscissa, as in Fig. and also in the 
ordinate. The departure from the quasi-Bloch theorem 
becomes very pronounced at the band edge. Above Vc, 
the participation ratios for states near the band edges 
decrease with increasing L, as shown in Fig. |^, but true 
exponential localization was not observed. There may 
be critical behavior, but more detailed analyses would 
be required to establish that. 

Participation ratios for part of the spectrum of the LJ 
chain are shown in Fig. |^. The behavior is very different 
from that for the FK chain, even though the behavior 
of the eigenvalue distribution of the two chains is simi- 
lar. There is a pronounced high-frequency band of the 
LJ chain where all eigenfunctions are exponentially lo- 
calized regardless of the smoothness of the underlying 
ground state, i.e., localization occurs for both V > Vc 



and V < Vc- The localized upper band is by no measure 
small: the band contains (Golden Mean)~^ ~ 38% of 
all states. The construction L x PR in Fig. ^ converges 
with successive Fibonacci approximations to distinct val- 
ues that might be called localization lengths. In fact, 
and in contrast to the results for the FK chain, examina- 
tion of the corresponding eigenfunctions shows that the 
probability distributions ipn very narrow and sharply 
peaked. They extend only a few atomic spacings and 
their participation ratios seem well-specified. The up- 
per branch in Fig. ^(b) is precisely a factor of two larger 
than the lower branch and corresponds to eigenfunctions 
that are concentrated at two spatial locations rather than 
one. What we term the upper band is itself a collection of 
smaller permitted and forbidden bands with gaps at suc- 
cessively finer scales as the Golden Mean is more closely 
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FIG. 6. Quasi-Bloch functions Uk{x) for the LJ chain 
for L = 377 at various quasimomenta k, with labelings as in 
Fig. |^(b). The values of k and the corresponding eigenenergies 
are (a) k = 0.38, E = 1.75 and 0.91 and (b) 0.42, 2.03 and 
0.93, for Vg = —80 and —110, respectively. These values 
are chosen from the "extended" region of Fig. ^ Above the 
transition, the functions are smeared out, as in Fig. ^ for the 
FK chain, again correlated with the participation ratio. 



with models which correspond to hard-core interactions, 
modest changes to compressible nearest-neighbor inter- 
actions may have drastic effects on phonon localization 
behaviour. 
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approximated. The lower band has more similarity to 
the FK case, both in the trends for the participation ra- 
tios and in the behavior of the quasi-Bloch functions, 
examples of which are shown in Fig. |6| However, there 
are some differences between the quasi-Bloch functions of 
the two models: the LJ quasi-Bloch functions shown for 
V < Vc have more oscillatory features than those of the 
FK functions, and the clumping along the abcissa and 
ordinate for V > Vc is so much more pronounced in the 
LJ chain than in the FK chain that the remnant of the 
underlying functions can barely be discerned. 

In conclusion, true, strong, phonon localization is ob- 
served in the upper band of the nearest-neighbor LJ 
chain. In fact the observation of a structure interpreted 
as a vibration within a domain wall for a uniaxially mod- 
ulated 2D lattice ||2^ was a motivation for examining 
the ID LJ chain. However, this localization is present 
regardless of the changes in the level statistics and pin- 
ning/unpinning of the chain, i.e., the transition by break- 
ing of analyticity. The behaviour of a closely related 
model, the FK chain, is strikingly different: the transi- 
tion by breaking of analyticity does affect localization in 
the system. These results are in accord with the widely 
held view that quasiperiodic chains lie on the border be- 
tween localization and delocalization and, thus, are very 
sensitive to even minor perturbations. Not only is it in- 
adequate to study localization in quasiperiodic systems 
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